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METHOD AND APPARATUS FOR 
TIME-FREE PROCESSING OF GPS SIGNALS 

BACKGROUND OF THE DISCLOSURE 

5 

1. Field of the Invention 

The invention relates to satellite-based position 
location and, more particularly, the invention relates to a 
method and apparatus for time-free processing of global 
10 positioning system (GPS ) signals. 

2 . Description of the Background Art 

Global Positioning System (GPS) receivers normally 
determine their position by computing time delays between 
15 transmission and reception of signals transmitted from 
satellites and received by the receiver on or near the 
surface of the earth. The time delays multiplied by the 
speed of light provides the distance from the receiver to 
each of the satellites that are in view of the receiver. 

2 0 The GPS satellites transmit to the receivers satellite- 

positioning data, so called "ephemeris" data. In addition 
to the ephemeris data, the satellites transmit to the 
receiver absolute time information associated with the 
satellite signal, i.e., the absolute time signal is sent as 
25 a second of the week signal. This absolute time signal 
allows the receiver to unambiguously determine a time tag 
for when each received signal was transmitted by each 
satellite. By knowing the exact time of transmission of 
each of the signals, the receiver uses the ephemeris data 

3 0 to calculate where each satellite was when it transmitted a 

signal. Finally, the receiver combines the knowledge of 
satellite positions with the computed distances to the 
satellites to compute the receiver position. 

The process of searching for and acquiring GPS 
35 signals, and reading the ephemeris and related data, 

including absolute time, for a multiplicity of satellites 
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is time consuming and introduces unacceptable delays m 
computing the receiver position. In addition, in many 
situations, there may be blockage of the satellite signals. 
In these cases the received signal level can be too low to 
demodulate and derive the satellite data without error. 
However, in these degraded signal situations, the receiver 
is capable of tracking the satellite signals, and measuring 
time delays (and hence distance) , if an external source of 
ephemeris and absolute time is available. 

Several innovations have been made to provide "GPS 
Aiding" that consists of external sources of ephemeris (or 
equivalent) data and absolute time information. The aiding 
information is transmitted to the GPS receiver using some 
alternative form of communication (usually wireless, such 
as cellular data channels) . Thanks to the use of GPS 
Aiding, GPS receivers can operate in areas where signal 
levels are too low for traditional GPS to function 
properly. 

All GPS Aiding, thus far invented, requires accurate 
external knowledge of the absolute time, so that the 
satellite positions can be accurately determined. The 
absolute time is required to an accuracy of between 1 
millisecond and 10 milliseconds. Unfortunately, there are 
desired implementations of GPS Aiding where absolute time 
cannot easily be obtained to this accuracy at the GPS 
receiver. For example: the AMPS cellular phone system does 
not support time information; nor (currently) does the 
North American TDMA cellular phone system. The GSM 
cellular phone system supports timing information, but may 
have different time references in different geographical 
regions. In these situations it is desirable to provide a 
method for computing GPS receiver position without knowing 
the absolute time. 

More specifically, Global Positioning System (GPS) 
receivers receive GPS signals transmitted from orbiting GPi 
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satellites containing unique pseudo-random noise (PN) 
codes. The GPS receivers determine the time delays between 
transmission and reception of the signals by comparing time 
shifts between the received PN code signal sequence and 
5 internally generated PN signal sequences. 

Each transmitted GPS signal is a direct sequence 
spread spectrum signal. The signals available for 
commercial use are provided by the Standard Positioning 
Service. These signals utilize a direct sequence spreading 
10 signal with a 1.023 MHz spread rate on a carrier at 1575.42 
MHz (the LI frequency) . Each satellite transmits a unique 
PN code (known as the C/A code) that identifies the 
particular satellite, and allows signals transmitted 
simultaneously from several satellites to be received 
15 simultaneously by a receiver with very little interference 
of any one signal by another. The PN code sequence length 
is 1023 chips, corresponding to a 1 millisecond time 
period. One cycle of 1023 chips is called a PN frame. Each 
received GPS signal is constructed from the 1.023 MHz 
20 repetitive PN pattern of 1023 chips. At very low signal 
levels the PN pattern may still be observed, to provide 
unambiguous time delay measurements, by processing, and 
essentially averaging, many PN frames. These measured time 
delays are called "sub-millisecond pseudoranges" , since 

2 5 they are known modulo the 1 millisecond PN frame 

boundaries. Once the absolute time delay can be 
calculated, by resolving the integer number of milliseconds 
associated with each delay to each satellite, then one has 
true, unambiguous, pseudoranges. The process of resolving 

3 0 the unambiguous pseudoranges is known as "integer 

millisecond ambiguity resolution" . 

A set of four pseudoranges together with a knowledge 
of the absolute times of transmissions of the GPS signals 
and satellite positions at those absolute times is 
35 sufficient to solve for the position of the GPS receiver. 
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The absolute times of transmission are broadcast from the 
satellites in the Navigation Message, which is superimposed 
on the 1.023 MHz PN code at a lower, 50Hz, data rate. This 
50 Hz signal is a binary phase shift keyed (BPSK) data 
stream with bit boundaries aligned with the beginning of 
the PN frame. There are exactly 20 PN frames per data bit 
period (20 milliseconds) . The 50 Hz signal contains data 
bits describing the GPS satellite orbits, satellite clock 
corrections, time of week information, and other system 
parameters . 

The absolute times associated with the satellite 
transmissions are determined in conventional GPS receivers 
by reading the Time of Week (TOW) data in the Navigation 
Message of the GPS signal. In the standard method of time 
determination, a conventional GPS receiver decodes and 
synchronizes to the 50 baud data stream. The 50 baud 
signal is arranged into 3 0 -bit words grouped into subframes 
of 10 words, with a length of 300 bits and a duration of 
six seconds. Five subframes comprise a frame of 1500 bits 
and a duration of 3 0 seconds, and 2 5 frames comprise a 
superframe with a duration of 12.5 minutes. A superframe 
contains the complete Navigation Message. The data bit 
subframes that occur every six seconds contain bits that 
provide the TOW to six second resolution. The 50 baud data 
stream is aligned with the C/A code transitions so that the 
arrival time of a data bit edge (on a 2 0 ms interval) 
resolves the absolute transmission time to the nearest 2 0 
milliseconds. Precision synchronization to bit boundaries 
can resolve the absolute transmission time to one 
millisecond or less. 

The absolute times associated with the satellite 
signals are determined in Wireless Aided-GPS receivers by 
having an external timing source that is calibrated to GPS 
time then using this time to provide a precise time tag at 
the time of reception of the satellite signal. The time of 
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reception minus the pseudorange gives the absolute time of 
transmission for each satellite (with the pseudorange 
expressed in time units, reflecting the transmission- 
reception time delay) . 

The absolute times of transmission are needed in order 
to determine the positions of the satellites at the times 
of transmission and hence to determine the position of the 
GPS receiver. GPS satellites move at approximately 3.9 
km/s, and thus the range of the satellite, observed from 
the earth, changes at a rate of at most ± 800 m/s. Absolute 
timing errors result in range errors of up to 0 . 8 m for 
each millisecond of timing error. These range errors 
produce a similarly sized error in the GPS receiver 
position. Hence, absolute time accuracy of 10 ms is 
sufficient for position accuracy of approximately 10m. 
Absolute timing errors of much more than 10ms will result 
in large position errors, and so typical GPS receivers have 
required absolute time to approximately 10 millisecond 
accuracy or better. 

Note that absolute timing errors also introduce errors 
as a result of the GPS satellite clock drift, but these are 
so much smaller than the satellite position error that they 
can be ignored for the purposes of this explanation (GPS 
clocks drift typically less than 0.1 nanoseconds per 
second, and the observed range to the satellite is affected 
by the GPS clock drift multiplied by the speed of light, 
this error is less than 0.03 m/s, about 25 thousand times 
smaller than errors caused by changes in satellite 
position) . 

There is another time parameter closely associated 
with GPS positioning, this is the sub-millisecond offset in 
the time reference used to measure the sub-millisecond 
pseudorange. This offset affects all the measurements 
equally, and for this reason it is known as the "common 
mode error" . 
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The common mode error should not be confused with the 
absolute time error. As discussed above, an absolute time 
error of 1 millisecond leads to range errors of up to 0.8 
meters while an absolute time error of 1 microsecond would 
cause an almost unobservable range error of less than 1 
millimeter. However, a common mode error of 1 microsecond 
causes a pseudorange error of 1 microsecond multiplied by 
the speed of light, that is 300 meters. 

Because common mode errors have such a large effect on 
pseudoranges, and because it is practically very difficult 
to calibrate the common mode error, traditional GPS 
receivers treat the common mode error as an unknown that 
must be solved for, along with position, once sufficiently 
many pseudoranges have been measured at a particular 
receiver. However, no traditional GPS receivers solve for 
absolute time error instead relying on the fact that they 
know absolute time to the required accuracy (of 10 
milliseconds or better) . 

Therefore, a need exists in the art for a method and 
apparatus that processes GPS signals without using absolute 
time . 

SUMMARY OF THE INVENTION 

The present invention is a method and apparatus for 
computing GPS receiver position without using absolute time 
information transmitted by a satellite or by an alternative 
source of timing available at the GPS receiver. In an 
embodiment of the invention, the GPS receiver is contained 
in an integrated receiver that also includes a wireless 
communication transceiver, but does not have access to an 
accurate source of absolute time information. The wireless 
transceiver communicates through a wireless network to a 
server. The GPS receiver measures satellite pseudoranges 
and uses the wireless communication transceiver to send the 
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pseudoranges to the server. The server fits the 
pseudoranges to a mathematical model in which the GPS 
receiver position and the absolute time are unknown 
parameters. The server then computes a position and 
absolute time that best fit the model, thus yielding the 
correct position for the GPS receiver, and the absolute 
time at which the pseudorange measurements were made. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The teachings of the present invention can be readily 
understood by considering the following detailed 
description in conjunction with the accompanying drawings, 
in which: 

FIG. 1 depicts a block diagram of apparatus for 
computing a GPS receiver location without knowledge of 
absolute time; 

FIG. 2 depicts a flow diagram representing the 
operation of the apparatus of FIG. 1; 

FIG. 3 depicts a flow diagram representing a method of 
computing pseudoranges in accordance with the invention; 

FIG. 4 depicts a flow diagram representing a method of 
computing receiver position and absolute time in an 
alternative embodiment of the invention; 

FIGS. 5A and 5B graphically depict a grid (5A) used 
for producing the residual magnitudes (5B) of position 
error; and 

FIG. 6 depicts a graph of the residual time error used 
in computing an absolute time. 

To facilitate understanding, identical reference 
numerals have been used, where possible, to designate 
identical elements that are common to the figures. 
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DETAILED DESCRIPTION 

The invention is a method and apparatus for 
determining position and time in a global positioning 
system (GPS) without having access, at the GPS receiver, to 
absolute time information. In the following description, 
for purposes of explanation, numerous specific details are 
set forth in order to provide a thorough understanding of 
the present invention. It will be evident, however, to one 
skilled in the art that the present invention may be 
practiced without these specific details. 

FIG. 1 depicts one embodiment of the present invention 
comprising an integrated mobile receiver 102 coupled to a 
server 121 via a wireless link 150. A GPS receiver 108 is 
contained in the integrated receiver 102 along with a 
wireless communication transceiver 112. The GPS receiver 
108 measures only sub-millisecond pseudoranges with respect 
to GPS satellites that are in view of the receiver 108, and 
then sends these sub-millisecond pseudo-ranges to a server 
121 using a wireless communication link 150. The server 
121 forms an approximate, a-priori estimate of the GPS 
receiver position from the known location of a wireless 
tower 118 used to receive the wireless communication. The 
server 121 also allocates a time tag from its own real time 
clock, thus creating an a-priori estimate of the absolute 
time at which the GPS receiver 10 8 received the GPS signals 
from the satellites. If the a-priori position is within 
100 km of the true position, and the a-priori absolute time 
of reception is within one minute of the true (unknown) 
time of reception, then the server 121 can resolve the 
integer milliseconds, thus turning the sub-millisecond 
pseudoranges into true pseudoranges. 

Next, the server 121 uses the pseudoranges to solve 
for the unknown receiver position and absolute time. The 
server comprises a central processing unit (CPU) 122, a 
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server clock 124, a tower location database 12 8, CPU 
support circuits 152, and memory 154. The support circuits 
comprise well-known circuits that facilitate operation of 
the CPU such as clock circuits, cache, power supplies, I/O 
interface circuits, and the like. The memory 154 may be 
random access memory, read only memory, removable storage, 
hard disk storage, or any combination of these memory 
devices . 

In one embodiment of the invention, the common mode 
error is assumed to be totally unknown at the server 121. 
In the one embodiment of the invention, the server 121 
assumes an a-priori common mode error of zero, although it 
will be understood that any other a-priori common mode 
error could be used, with no change in the results. With 
the five a-priori estimates of the unknown parameters (3 
coordinates of position, 1 absolute time, 1 common mode 
error) the server 121 creates a mathematical model relating 
the measured pseudoranges and a-priori information to the 
unknown parameters. The mathematical model can be written 
as a linear equation, which, when solved, yields the 
correct position and absolute time. 

More specifically, GPS signals 104 from a plurality of 
satellites (not shown) are received at the GPS antenna 106. 
The received signals are coupled to the GPS receiver 108. 
The GPS receiver 108 processes the GPS signals to form sub- 
millisecond pseudoranges on path 110, which are coupled to 
the communication transceiver 112 and transmitted from the 
communication antenna 116 through a wireless network such 
as a cellular telephone network. The transmission from the 
integrated receiver 102 is received by a nearby radio tower 
118, e.g., a cellular telephone tower. The sub-millisecond 
pseudoranges and the radio tower ID are sent from said 
radio tower 118 to the server 121. In the server 121 the 
server clock 124 is used to provide a time-tag when the 
sub-millisecond pseudoranges are received at the server. 
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The server 121 passes the radio tower ID along path 12 6 to 
a tower location database 12 8, and extracts a location for 
the tower from the database 128. The tower location is 
coupled to the CPU 122 along path 130. 

The satellite ephemeris data, for all satellites 
represented by the sub-millisecond pseudorange data, is 
provided to the server from some external source 12 5 (such 
as another GPS receiver located in the vicinity of the 
server with a clear view of the sky, or some other source 
such as a network of GPS receivers) . Note that, for 
simplicity, the term "ephemeris" is used to mean the 
satellite orbital parameters, as well as the satellite 
clock parameters. The CPU 122 of the server 121 combines 
the sub-millisecond pseudoranges , radio tower location, 
server time, and ephemeris to form the correct GPS receiver 
position and absolute time of reception of signals at the 
GPS receiver 108. 

The apparatus above assumes that the GPS receiver 108 
is not capable of reliably receiving the absolute time 
information and ephemeris data, i.e., the GPS receiver is 
used indoors, such that the processing of ephemeris data 
is accomplished in the server 121. However, in some 
instances, rather than having the pseudo range data 
supplied to the server, the server (or some other source) 
can supply the ephemeris data and clock signal to the 
mobile device 102 and the mobile device can perform the 
position calculation. In such an embodiment of the 
invention, a CPU and clock (similar to 122 and 124) are 
located in the mobile device 102 to facilitate signal 
processing in the same manner as is described with respect 
to the server 121. 

FIG. 2 is a flowchart illustration of the process 200 
that is performed by the server CPU 122 of FIG. 1. At step 
2 02, the server clock signal is used to provide an a-priori 
estimate of the absolute time of reception of the GPS 
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signals at the GPS receiver. It will be understood that 
the use of the server clock is one embodiment used to 
exemplify the current invention and, in general, the a- 
priori estimate of time could come from a time source other 
than the server clock. The present invention is applicable 
irrespective of the source, or quality, of the a-priori 
estimate of absolute time. To simplify the exposition, this 
particular embodiment is assumed to have a server clock 
that provides a time tag within one minute of the actual 
absolute time of reception of the GPS signals at the GPS 
receiver. It will be understood that this . simplifying 
assumption, while often true in practice, is not a 
necessary part of the invention, and has been used here 
only to simplify the explanation of the invention. Later 
in this specification, this simplifying assumption is 
removed . 

At step 206, the tower location is provided to the CPU 
as an a-priori estimate of the GPS receiver position. It 
will be understood that the use of the tower location is 
just one embodiment of any number of a-priori positions 
that could be used (for example: a previously calculated 
position for the same GPS receiver 108 could be used as an 
a-priori position, or a combination of positions of 
recently used towers, or the a-priori position could simply 
be guessed) . The present invention is applicable 
irrespective of the source, or quality, of that a-priori 
position. To simplify the exposition, this particular 
embodiment is assumed to have an a-priori position that is 
within 100 km of the true position of the GPS receiver 108. 
It will be understood that this simplifying assumption, 
while often true in practice, is not a necessary part of 
the invention, and has been used here only to simplify the 
explanation of the invention. Later in this specification 
this simplifying assumption is removed. 
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At step 2 04 and 2 08, the sub-millisecond pseudoranges 
and ephemeris for the appropriate satellites that are in 
view of GPS receiver are also provided as inputs to the 
process 200. 

At step 210, the sub-millisecond pseudorange integers 
are resolved by a process described below with respect to 
FIG. 3. Having resolved the sub-millisecond pseudorange 
integers, the process creates full pseudoranges. 

At step 212, the expected pseudoranges are formed. 
These expected pseudoranges are the pseudoranges that woul 
be measured if all the a-priori parameters (a-priori 
position, a-priori absolute time of measurement, and a- 
priori common mode error) were in fact the actual values o 
these parameters . The expected pseudoranges are denoted 
t ± , the index i denotes the appropriate satellite. 

At step 214, the a-priori pseudorange residuals are 
formed, these residuals are defined as the difference 
between the measured pseudoranges (denoted p, ) and the 
expected pseudoranges (rj . The a-priori pseudorange 
residuals are denoted u i . 

At step 216 a mathematical model is formed, relating 
to x, where u is a vector of u £ and x is a vector of the 
updates to the a-priori values of position, common-mode 
error, and absolute time of reception: 



where n is the number of pseudoranges. The pseudoranges 
are expressed in units of length (e.g. meters) . 



-13- 



The position updates, x, y, z, are expressed in units of 
length (e.g. meters) and the time updates t c , t s are in 
units of time (e.g. seconds) 

An embodiment of a mathematical model relating these 
5 two vectors is a Taylor series, where the first term in the 
series is the first derivative of u with respect to x, the 
second term contains the second derivative, and so on. In 
one embodiment of the process, the invention uses a 
linearized model that keeps only the first derivative in 
10 the Taylor series. This gives the following equation 
relating u to x: 



dp x ldx dpjdy dpjdz dpjdtc dp x ldt s 
dpjdx dpjdy dpjdz dpjdt c dpjdt s 



dpjdx dpjdy dpjdz c - p, 
dp n /dx dp n ldy dp n /dz c -p n \ 



The particular values of dpjdx , dpjdy , and dpjdz depend 
15 on the coordinate system used to describe the a-priori 
position. These terms in the first three columns of the 
matrix H are well known in the art, and further explanation 
is not required. The fourth column of the matrix is c, the 
speed of light, this part of the model is also standard in 
20 the art. The novel aspect of this invention requires the 
inclusion of the fifth column in the matrix. This fifth 
column exactly models the relationship between the unknown 
error in the a-priori estimate of absolute time, and the 
measured pseudoranges . Furthermore the terms in this 
2 5 column are the rate of change of the pseudoranges with 
respect to time and can be exactly calculated from the 
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ephemeris data. Hence every term in the matrix H is known 
and, provided there are five or more pseudoranges available 
at the GPS receiver, the values of x can be calculated 
using linear algebra. 

At step 22 0, the GPS receiver position is computed by 
adding the updates x,y,z, to the a-priori position, and the 
absolute time of reception is formed by adding the update 
t s to the a-priori time of reception. If the a-priori 
position and a-priori absolute time were close enough to 
the true position and true absolute time, then one pass 
through the process 200 will yield results to the required 
accuracy. However, if the first pass through the process 
2 00 does not immediately converge on the required accuracy, 
then the result 222 is used to form a new a-priori time of 
reception estimate for step 2 02 and a new a-priori 
position estimate for step 206, and the process 200 is 
repeated until the result converges on the correct result 
(typically very few iterations are required because the 
linearization using the first order Taylor series is a very 
good mathematical description of the complete, non-linear, 
system, thanks to the fact that the satellite ranges are so 
much further from the earth than the error in the a-priori 
position) . 

It will be understood that the Taylor series is just one 
example of a mathematical model relating the unknown 
position and absolute time to the measured pseudoranges. 
The present invention is equally valid with other models, 
such as non- linear models, which can be solved through 
techniques such as iteratively fitting the unknown 
parameters until an acceptable solution is obtained. 

If, as assumed above, the a-priori position and 
a-priori absolute time are within 100km and 1 minute 
respectively, then the result 222 will be correct. 
However, if the a-priori position and time are not known 
within these limits, then the incorrect integers may be 
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formed at step 210 and the incorrect result 222 may be 
obtained. In this case the a-posteriori residuals, formed 
at step 224, will be used, as described below with respect 
to FIG. 4, to detect this errant condition and then 
B different a-priori values will be used. 

FIG. 3 is a flowchart of an illustrative process 300 
that resolves the sub-millisecond pseudorange integers 
(step 210 of FIG. 2) . To simplify the exposition, this 
particular embodiment is assumed to have an a-priori 

10 position that is within 100 km of the true position of the 
GPS receiver, and an a-priori absolute time estimate that 
is within one minute of the true absolute time of reception 
at the GPS receiver. It will be understood that these 
simplifying assumptions, while often true in practice, are 

15 not a necessary part of the invention, and have been used 
here only to simplify the explanation of this embodiment of 
the invention. In the description relating to FIG. 4, FIG. 
5, and FIG. 6, these simplifying assumptions are removed. 
At step 3 08, the process 3 00 calculates the expected 

2 0 pseudoranges using ephemeris data (provided in step 2 08) 

for the satellites along with the a-priori absolute time of 
reception (provided in step 2 02) and the a-priori position 
(provided in step 206) . As done throughout this 
specification, the term ephemeris is used to mean the 
25 satellite orbital parameters as well as the satellite clock 
parameters . 

At step 310, a single satellite is chosen as a 
reference satellite. In the preferred embodiment the 
satellite with the highest elevation angle (from the a- 

3 0 priori position) is chosen as the reference, but it will be 

understood that it is not important which satellite is used 
as the reference. The expected pseudorange for the 
reference satellite is denoted r 0 (path 312). The expected 
pseudoranges for the other satellites are denoted r x (path 
35 314) . 
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At step 318, an integer is assigned to the reference 
satellite. The integer must satisfy the equation: 



where c is the speed of light, expressed in m/s, t c is the 
common mode error, and e 0 is the error in the expected 
pseudorange introduced by the combined errors in the a- 
priori position and a-priori absolute time. 

Those skilled in the art will understand that an 
arbitrary integer may be assigned, since in the following 
discussion it is seen that the common mode error will 
absorb any errors made in this integer, as long as exactly 
equal errors are made in all the other integers. The 
integer N 0 is assigned according to the following equation: 



At step 322, the integer millisecond values for the 
remaining satellites are calculated using the sub- 
millisecond pseudoranges 320 for all the satellites along 
with N 0 . Whatever value of N 0 was chosen above implies an 
associated common mode error t c . The values of N ± are 
chosen to satisfy the following equation, which relates N if 
the measured sub-millisecond pseudorange (sj , the expected 
pseudorange (r ± ) , and the common mode error (t c ) . 



where, e ± is the error in the expected pseudorange 
introduced by the combined errors in the a-priori position 
and a-priori absolute time. In the preferred approach, the 
corresponding equation for N 0 is subtracted from the above 
equation, this exactly cancels the term t c , since this term 
is common (by definition) for all satellites. This yields 
the following expression for N t 



The above equations represent one embodiment of a process 
to calculate the integers. It will be understood that any 
expression may be used to compute these integers, provided 



N 0 *c/10 3 + 5 0 - t c =r 0 -e 0 , 



N Q = round ((r 0 - s Q )* 1 0 3 / c) . 



N, * c/10 3 + s, -t c =r,-e, r 
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the relationship to t c is consistently maintained for all 
the integers . 

In the above description the assumption was made that 
the a-priori position was within 100km of the true 
5 position, and the a-priori absolute time was within 1 
minute of the true time. For all GPS satellites, the 
maximum pseudorange rate is ±800 m/s. Thus the maximum 
value of the error term e i will be 100km + 60s*0.8km/s = 
148 km. This is less than half of one C/A code millisecond 

10 epoch (i.e., less than half of one integer millisecond) and 
the rounding operation used above will always yield the 
correct integers. In the disclosure with respect to FIG. 
4, FIG. 5, FIG. 6, these two restrictions on a-priori 
position accuracy and a-priori time accuracy are removed. 

15 If the a-priori position is not known to within 100km, 

it will nonetheless be known to some finite accuracy limit. 
Similarly, if the a-priori absolute time is not known to 
within 1 minute it will be nonetheless be known to some 
finite accuracy limit. As described in the processes 200 

20 and 300, any a-priori position and time estimate within the 
100km and 1 minute constraints will yield the correct 
integers, correct GPS receiver position, and correct 
absolute time. In an embodiment of the current invention, 
the space of all possible a-priori positions is divided 

25 into a lOOkmXlOOkm lat-lon grid, with altitude assigned 

from a look-up table of topographical altitudes. Similarly 
the space of all possible a-priori absolute times is 
divided into 1 -minute segments. This yields a set of all 
possible a-priori positions and times. The process 200 is 

30 used iteratively, with each of the possible a-priori 

position and time values from the set. When an a-priori 
position and time is found within 100km and 1 minute of 
truth, then the a-posteriori residuals will be small, and 
the correct GPS receiver position and absolute time will be 

35 computed, as described above. 
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An embodiment of this process 400 is shown in FIG. 4. 
At step 402, all possible a-priori positions and residuals 
are formed into a set. In one embodiment the set is 
organized as a lOOkmXlOOkm grid, with altitudes assigned 
5 from a lookup table of topographical heights, and with time 
segmented into 1 -minute segments. This is a convenient 
representation because, as discussed above, the combination 
of a 100km grid plus the maximum pseudorange rate times 1- 
minute gives a maximum possible estimated error of just 

10 less than half a millisecond of range, which is necessary 
for the integer resolution process 300 to select the 
correct integers. However it will be understood that any 
number of different methods may be used to organize the set 
of all possible a-priori values, including creating new 

15 elements of the set dynamically, based on the results 
obtained from previously used values. 

At step 404, the process 400 selects one possible a- 
priori position and time combination. These values are 
used in the process 200. At step 406, the process 400 

2 0 examines the a-posteriori residuals that are produced in 
step 224 of FIG. 2. If the correct GPS receiver position 
and absolute time have been calculated, then the magnitude 
of the residuals will be small (that is, of the same order 
as the pseudorange measurement errors - tens of meters) . 

2 5 If the a-priori position and time were far enough away from 

the truth that the integer ambiguities were not correctly 
resolved, then the residuals will be large (that is, of the 
order of one millisecond epoch - many kilometers) . If the 
residuals are large then the candidate a-priori position 

3 0 and time are incorrect and they are removed from the set of 

possibilities. The process 400 is iterated until the 
correct position and absolute time are calculated. 

It will be understood by those skilled in the art that 
even if the absolute time is available at the GPS receiver, 
35 the process of integer ambiguity resolution has 
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traditionally required an initial estimate of position 
close enough to the true position for the integers to be 
uniquely defined by said initial estimate. The current 
invention provides a novel means of computing the correct 
integers without requiring an accurate initial estimate of 
position . 

FIGS. 5A and 5B respectively depict a grid 502 used to 
determine receiver position in an embodiment of the 
invention and the residual magnitudes 5 06 computed during 
the position calculation process 200. In this example, the 
a-priori position is assigned as an arbitrary guess, in the 
middle of North America. Then the process 2 00 is performed 
for each possible a-priori position (grid point 504) , and 
the magnitude of the a-posteriori residuals is examined. 
As each wrong candidate is rejected another candidate 
(another grid point 504) is generated by searching outwards 
on a 1 degree x 1 degree grid 5 02. (Note that this 1 degree 
x 1 degree grid is a slightly different embodiment than the 
100km X 100km grid described earlier; both embodiments are 
guaranteed to give at least one a-priori position that will 
yield the correct integers, and hence the correct position 
and absolute time.) The a-priori position is completed by 
assigning an a-priori altitude from a lookup table of 
topographical heights. FIG 5A shows the 1,663 wrong 
candidates on a grid 504 and FIG 5B shows the corresponding 
residual magnitudes 506, each one of the order of an 
incorrect millisecond integer (i.e., many kilometers). 
Once the search reaches an a-priori position in the 
vicinity of the true position (in San Jose, CA) the 
mathematical model causes the correct result to "snap" into 
place, and the correct position and time are calculated. 
The a-priori position candidate number 1,664 (residual 
magnitude 508 and grid point 510) is approximately 175 km 
east of the true position, which, in this example is close 
enough for the position and time solution to "snap" into 
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place. The correct solution yields a residual that is 
approximately 3 0 meters, which is from one thousand to ten 
thousand times smaller than the incorrect residuals. 

The large difference between "small" residuals (tens 
5 of meters) and "large" residuals (tens to hundreds of 
kilometers) makes this embodiment work very well in 
practice. However, it will be understood that other 
methods may be used to test the quality of the result, 
including comparing the calculated position and absolute 

10 time to a position and time obtained through some other 
means - such as the location of the radio tower used in a 
wireless-aided system. It will also be appreciated that in 
order to have non-zero residuals, it is necessary to have 
more independent observations than unknowns . In the 

15 embodiments described thus far, there have been five 

unknown parameters: three coordinates of position, common 
mode error, and absolute time. Thus at least six 
independent measurements are required to have non-zero 
residuals. If six GPS satellites can be observed, then the 

2 0 six measurements can be obtained from them. If there are 
not six satellites observable then there are several steps 
that can be taken, many of which are standard in the art. 
The number of measurements may be increased by including 
measurements from other sources (such as range measurements 

2 5 based on time-of -arrival from a wireless system, or angle- 

of-arrival measurements measured in a wireless system, or 
any other independent measurement that can be obtained) . 

The number of observables can also be increased by 
including, as "pseudo-measurements", known constraints on 

3 0 position, for example known or approximately known altitude 

can be introduced into the mathematical model as a pseudo- 
measurement. In the embodiment specified above, where the 
mathematical model is represented by the equation u=i/x / a 
pseudo-measurement for known altitude may be created by 
35 first specifying the a-priori position in coordinates of 
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latitude, longitude, and altitude, then by setting the a- 
priori altitude to the known altitude, then by adding a new- 
line to the matrix equation: 



5 This approach effectively adds another measurement or 
observable to the mathematical model. This approach is 
standard in the art, and it is understood that it applies 
to any constraints that may be useful in solving for the 
unknown parameters . 

10 Another approach is to reduce the number of unknown 

parameters. This may be performed by removing known, or 
approximately known parameters. The most commonly known 
parameter is altitude, and it can be removed from the 
mathematical model. Similarly the common mode error may be 

15 calibrated (for example, if the invention is implemented in 
a system with access to a stable oscillator) and removed 
from the mathematical model. 

It will be appreciated that many combinations of the 
disclosed techniques may be applied to compute unknown 

2 0 values, including unknown absolute time. 

For example, using the techniques of this invention, 
one can compute the time parameters alone, without 
computing the position. This is done, in the preferred 
embodiment, by fixing the position in the mathematical 

2 5 model to the a-priori position, and computing the remaining 
two unknown parameters : common mode error and absolute 
time . 

FIG 6. is a graphical representation 600 of the 
residuals magnitudes (axis 602) associated with the 
30 different a-priori times (axis 604) that were attempted in 
an example embodiment of the invention. In this particular 
example, a range of possible times, each one-minute apart, 
is attempted for each of the grid points shown in FIG. 5. 
The first a-priori absolute time was chosen by guessing a 
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time that turns out to be approximately two -and- a -ha If 
hours later than the true absolute time of reception. When 
the process 400 applies an a-priori position approximately 
175 km away from the true position, and an a-priori time 
within one minute of the true absolute time, the 
mathematical model calculates the correct position and 
time, as shown in FIG 5. In general, the mathematical model 
will calculate the correct position and time as soon as the 
a-priori position and time are close enough to cause 
process 300 to calculate the correct integers. As 
discussed, the preferred embodiment is guaranteed to find 
at least one such a-priori position and time, by creating 
an appropriate grid, and appropriately spaced time 
intervals . 

Although the methods and apparatus of the present 
invention have been described with reference to GPS 
satellites, it will be appreciated that the teachings are 
equally applicable to positioning systems which utilize 
pseudolites or a combination of satellites and pseudolites. 
Pseudolites are ground based transmitters that broadcast a 
PN code {similar to the GPS signal) which may be modulated 
on an L-band carrier signal, generally synchronized with 
GPS time, the term "satellite", as used herein, is intended 
to include pseudolites or equivalents of pseudolites, and 
the term "GPS signals", as used herein, is intended to 
include GPS-like signals from pseudolites or equivalents of 
pseudolites . 

In the preceding discussion the invention has been 
described with reference to application upon the United 
States Global Positioning System (GPS) . It should be 
evident, however, that these methods are equally applicable 
to similar satellite systems, and in particular, the 
Russian Glonass system and the European Galileo system. The 
term "GPS" used herein includes such alternative satellite 
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positioning systems, including the Russian Glonass system 
and the European Galileo system. 

Although various embodiments which incorporate the 
teachings of the present invention have been shown and 
described in detail herein, those skilled in the art can 
readily devise many other varied embodiments that still 
incorporate these teachings. 
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What is claimed is: 

5 1. A method for calculating position and time of a GPS 
receiver comprising: 

providing pseudoranges that estimate the range of the 
GPS receiver to a plurality of GPS satellites ; 

providing an estimate of an absolute time of reception 
10 of a plurality of satellite signals; 

providing an estimate of a position of the GPS 
receiver ; 

providing satellite ephemeris data; 

computing absolute position and absolute time using 
15 said pseudoranges by updating said estimate of an absolute 
time and the estimate of position of the GPS receiver. 

2 . The method of claim 1 wherein said pseudoreanges are 
sub-millisecond pseudoranges. 

20 

3. The method of claim 1, wherein a first instance of said 
estimate of absolute time is in error by more than 10 
milliseconds . 

25 4. The method of claim 1, wherein said estimate of absolute 
time is provided by a clock that is not linked to a GPS 
reference time. 

5. The method of claim 1, wherein said estimates of 
30 position and absolute time may be arbitrary guesses. 

6. The method of claim 1, further comprising: 

forming a-priori pseudorange-residuals ; and 
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forming a mathematical model relating said a-priori 
pseudorange-residuals to said updates of said position and 
absolute time estimates; and 

computing updates of the position and absolute time 
5 that fit the mathematical model. 

7. The method of claim 6, wherein said a-priori range- 
residuals are a difference between expected ranges, from 
said satellites to said a-priori position estimate, and the 

10 pseudoranges . 

8. The method of claim 6, wherein said expected ranges are 
computed at a time given by said a priori time estimate. 

15 9. The method of claim 6, wherein said mathematical model 
is a linearization of a Taylor series of a non-linear 
mathematical model. 

10. The method of claim 9, wherein said linearization is of 
2 0 the form: 

u, = \dp,/dx, dp, Idy, dpjdz, dp,/dt c , dp, / dt s ~\*h, y, z , t c , *sT 
where : 

2 5 Ui is the a-priori range-residual, for one satellite ; 

p z is the pseudorange for a satellite i; 

x, y and z are three coordinates of position updates; 

t c is a common mode error update; 

t s is an absolute time of reception update; and 

3 0 3 denotes the partial derivative. 

11. The method of claim 4, wherein one or more of the 
updates is assumed known, so that the remaining updates may 
be computed. 

35 
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12. The method of claim 10, wherein one or more of the 
updates x, y, z, or t c is assumed known, and set equal to 
an assumed value in the mathematical model, so that the 
remaining updates may be computed. 

5 

13. The method of claim 10, wherein one or more of the 
updates x, y, z, or t c is assumed known, and added to the 
model as a pseudo-measurement so that the remaining updates 
may be computed. 

10 

14. A process as in claim 6, wherein other measurements or 
constraints are used in the mathematical model. 

15. A process as in claim 6, wherein said a-priori 

15 position is obtained from a location of a radio tower used 
to communicate with a mobile device containing said GPS 
receiver . 

16. A process as in claim 6, wherein said absolute time 
20 estimate is obtained from a real time clock in a server, 

said server being located remotely from said GPS receiver. 

17 . A method for calculating absolute time of a GPS 
receiver comprising: 

2 5 providing pseudoranges that estimate the range of the 

GPS receiver to a plurality of GPS satellites; 

providing an estimate of position of the GPS receiver; 

and 

computing absolute time using the pseudoranges and the 

3 0 position estimate. 



18. The method of claim 1 wherein said pseudoranges are 
sub-millisecond pseudoranges. 
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19. The method of claim 17, further comprising: 
forming a-priori pseudorange-residuals ; and 
forming a mathematical model relating said a-priori 
5 pseudorange-residuals to updates of said estimate of 

position and absolute time; and 

computing updates of the position and absolute time 

that fit the mathematical model. 

10 20. The method of claim 19, wherein said a-priori range- 
residuals are a difference between expected ranges, from 
said satellites to said a-priori position estimate, and the 
pseudoranges . 

15 21. The method of claim 19, wherein said expected ranges 
are computed at a time given by said a priori time 
estimate . 

22. The method of claim 19, wherein said mathematical model 
20 is a linearization of a Taylor series of a non-linear 

mathematical model . 

23. The method of claim 22, wherein said linearization is 
of the form: 

25 

u, = \dp, fdx, dpjdy, dpjdz, dp,Idt c , dp,Jdt s ]*\x, y, z, t c , t s J 
where : 

u ± is the a-priori range-residual, for one satellite; 
30 p ± is the pseudorange for a satellite i; 

x, y and z are three coordinates of position updates ; 

t c is a common mode error update; 

t s is an absolute time of reception update; and 

d denotes the partial derivative. 

35 
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24. The method of claim 23, wherein one or more of the 
updates x, y, z, or t c is assumed known, and set equal to 
an assumed value in the mathematical model, so that the 
remaining updates may be computed. 

5 

25. The method of claim 24, wherein one or more of the 
updates x, y, z, or t c is assumed known, and added to the 
model as a pseudo-measurement so that the remaining updates 
may be computed. 

10 

26. A process as in claim 24, wherein other measurements 
or constraints are used in the mathematical model. 

27. A process as in claim 24, wherein said a-priori 

15 position is obtained from a location of a radio tower used 
to communicate with a mobile device containing said GPS 
receiver . 

28. A method of calculating a GPS position for a GPS 

20 receiver from partial pseudoranges that have ambiguity in a 
number of integer milliseconds, comprising: 

a) choosing an a-priori position of the GPS receiver; 

b) calculating integers conforming to said a-priori 
position; 

25 c) calculating a navigation solution; 

d) calculating a-posteriori residuals; and 

e) using a relative size of said a-posteriori 
residuals to determine if said calculated integers are 
correct ; and 

3 0 f ) repeating steps c) , d) and e) using another a- 

priori position until residuals having a magnitude below a 
predefined threshold are computed. 
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29. A method as in claim 28, wherein said a-priori 
position is not within a specified distance of the actual 
receiver position. 

5 30. A method as in claim 28, wherein said a-priori 

position is greater than 100km from the actual receiver 
position . 

31. A method as in claim 28, wherein said a-priori 

10 position is greater than one integer millisecond from the 
actual receiver position. 

32. A method as in claim 28, wherein said a-priori position 
is an arbitrary guess. 

15 

33 . A system for computing position and time for a GPS 
receiver comprising: 

a mobile device comprising a GPS receiver and a 
wireless transceiver; 
2 0 a server being in wireless communication with said 

mobile device; 

where said GPS receiver computes pseudoranges that estimate 
the range of the GPS receiver to a plurality of GPS 
satellites and the wireless transceiver transmits said 
2 5 pseudoranges to the server; 

where the server computes an absolute position and 
absolute time for the GPS receiver using the pseudoranges 
and an estimate of position and time. 

30 34. The system of claim 33 wherein said pseudoranges are 
sub-millisecond pseudoranges. 

35. The system of claim 33 wherein the position estimate 
is a position of a radio tower coupled to the server that 
35 receives signals from the wireless transceiver. 
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Abstract of the Disclosure 

5 A method and apparatus for computing GPS receiver 

position without using absolute time information 
transmitted by the satellite or by an alternative source of 
timing available at the GPS receiver. The GPS receiver is 
contained in an integrated receiver that also includes a 

10 wireless communication transceiver, but does not have 

access to an accurate source of absolute time information. 
The wireless transceiver communicates through a wireless 
network to a server. The GPS receiver measures satellite 
pseudoranges and uses the wireless communication 

15 transceiver to send the pseudoranges to the server. The 
server fits the pseudoranges to a mathematical model in 
which the GPS receiver position and the absolute time are 
unknown parameters. The server then computes a position 
and absolute time that best fit the model, thus yielding 

2 0 the correct position for the GPS receiver, and the absolute 
time at which the pseudorange measurements were made. 
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